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SURFACES USING VORTEX-LATTICE METHODOLOGY 
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NASA Langley Research Center 
SUMMARY 

A new subsonic method has been developed by which the mean camber surface 
can be determined for trimmed noncoplanar planforms with minimum vortex drag. 

This method uses a vortex lattice and overcomes previous difficulties with 
chord loading specification. This method uses a Trefftz plane analysis to 
determine the optimum span loading for minimum drag, then solves for the mean 
camber surface of the wing, which will provide the required loading. Pitching- 
moment or root-bending-moment constraints can be employed as well at the design 
lift coefficient. 

Sensitivity studies of vortex-lattice arrangement have been made with this 
method and are presented. Comparisons with other theories show generally good 
agreement. The versatility of the method is demonstrated by applying it to 
(1) isolated wings, (2) wing-canard configurations, (3) a tandem wing, and 
(4) a wing-winglet configuration. 

INTRODUCTION 

Configuration design for subsonic transports usually begins with the wing, 
after which the body and its effects are taken into account, and then the tails 
are sized and located by taking into account stability and control requirements. 

With the advent of highly maneuverable aircraft having closely coupled lifting 
surfaces, there has been an increased interest in changing the design order so 
that multiple surfaces could be designed together to yield a trimmed configura- 
tion with minimum induced drag at some specified lift coefficient. Such a com- 
bined design approach requires that the mutual interference of the lifting sur- 
faces be considered initially. 

Single planform design methods are available to optimize the mean camber 
surface, better called the local elevation surface, for wings flying at sub- 
sonic speeds (for example, ref. 1) and at supersonic speeds (for example, 
refs. 2 and 3). The design method presented in reference 1 was developed from 
an established analysis method (Multhopp type), also presented in reference 1, 
by using the same mathematical model, but the design method solves for the 
local mean slopes rather than the lifting pressures. In the usual implementa- 
tion of reference 1, the design lifting pressures are taken to be linear chord- 
wise, but must be represented in this solution by a sine series which oscillates 
about them. An example presented herein demonstrates that corresponding oscil- 
lations may appear in pressure distributions measured on wings which have been 
designed by the method of reference 1. The method developed herein overcomes 
this oscillatory lifting pressure behavior by specifying linear chord loadings 
at the outset. 
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The development approach used in the two-planform design problem will be 
similar to that used for a single planform. The analytic method employed, 
selected because of its geometric versatility, is the noncoplanar two-planform 
vortex-lattice method of reference 4. 

The design procedure is essentially an optimization or extremization prob- 
lem. Subsonic methods (for example, see refs. 5 and 6) are available for deter- 
mining the span load distributions on bent lifting lines in the Trefftz plane, 
but they do not describe the necessary local elevation surface. This is one of 
the objectives of the present method which will utilize the Lagrange multiplier 
technique (also employed in refs . 2 and 3) . The method of reference 4 is used 
to provide the needed geometrical relationships between the circulation and 
induced normal flow for complex planforms, as well as to compute the lift, drag, 
and pitching moment. 

This paper presents limited results of precision studies and comparisons 
with other methods and data and is a condensed version of reference 7. Several 
examples of solutions for configurations of recent interest are also presented. 

SYMBOLS 

A element of aerodynamic influence function matrix A which con- 

tains induced normal flow at 1th point due to nth horseshoe 

N N 

vortex of unit strength; total number of elements is y * — 
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aspect ratio 

fractional chord location where chord load changes from constant 
value to linearly varying value toward zero at trailing edge 

coefficients in spanwise scaling polynomial 

wing span 

drag coefficient 

drag coefficient at ■ 0 


vortex or induced drag coefficient. 
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X,Y,Z 

X,Y,Z 
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Ax 


section lift coefficient 
reference chord 

= p c a + 0.7 5^| (brackets indicate "take the greatest integer") 

maximum number of spanwise scaling terms in solution technique 
for wings without dihedral 

lift 

pitching moment about coordinate origin 
free-stream Mach number 

number of span stations where pressure modes are defined as used 
in reference 1 

maximum number of elemental panels on both sides of configura- 
tion; maximum number of chordal control points at each of m 
span stations as used in reference 1 

number of elemental panels from leading to trailing edge in 
chordwise row 

total number of (chordwise) rows in spanwise direction of 
elemental panels on configuration semispan 

free-stream dynamic pressure 
reference area 

horseshoe vortex semiwidth in plane of horseshoe (see fig. 1) 
free-stream velocity 

axis system of given horseshoe vortex (see fig. 1) 

body-axis system for planform (see fig. 1) 

wind-axis system for planform (see sketch (a)) 

distance along X-, Y-, and Z-axis, respectively 

distance along X-, Y-, and Z-axis, respectively 

incremental movement of X-Y coordinate origin in streamwise 
direction 





(3z/3x) 


y and z distances from image vortices located on right half 
of plane of symmetry, as viewed from behind, to points on ( 
left panel 

canard height with respect to wing plane, positive down 

local elevation normalized by local chord, referenced to local 
trailing-edge height, positive down 

Ith elemental local slope in vector {3z/3x} of N/2 elements 
(see eq. (1)) 

angle of attack, deg 

vortex strength of nth element in vector { r } of N/2 elements 

incidence angle, positive leading edge up, deg 

independent variable in extremization process 

nondimensional spanwise coordinate based on local planform 
semispan 

distance along local chord normalized by local chord 

fractional chordwise location of point where mean camber height 
is to be computed (see eq. (14)) 

dihedral angle from trailing vortex to point on left panel being 
influenced; a measured from left panel, o' measured from 
right panel 

horseshoe vortex dihedral angle in Y-Z plane on left wing panel, 
deg 

horseshoe vortex dihedral angle on right wing panel, 
deg 


Subscripts: 


canard 


design 

indices to vary over the range indicated 
leading edge 

associated with slope point and horseshoe vortex, respectively, 
ranging from 1 to N/2 

left trailing leg 
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R 


right trailing leg 
reference value 


ref 

w wing 

Matrix notation: 

{ } column vector 

[ } square matrix 




Flow angle of attack determined 
at each sic; e point 

Typical spanwise 
vortex filament 


U 


Vortex-lattice trailing filaments 

Sketch (a) 



THEORETICAL DEVELOPMENT 

This section presents the application of vortex-lattice methodology to the 
mean-camber-surface design of two lifting planforms which may be separated ver- 
tically and have dihedral. For a given planform, local vertical displacements 
of the surfaces with respect to their chord lines in the wing axis (see 
sketch (a)) are assumed to be negligible; however, vertical displacements of 
the solution surfaces due to planform separation or dihedral are included. The 
wakes of these bt.it lifting planforms are assumed to lie in their respective 
extended bent chord planes with no roll up. For a two-planform configuration 
the resulting local elevation surface solutions are those for which both the 
vortex drag is minimized at the design lift coefficient and the pitching moment 
is constrained to be zero about the origin. For an isolated planform no 
pitching-moment constraint is imposed. Thus, the solution is the local eleva- 
tion surface yielding the minimum vortex drag at the design lift coefficient. 
Lagrange multipliers together with suitable interpolating and integrating pro- 
cedures are used to obtain the solutions. The details of the solution are given 
in the following five subsections. 


Relationship Between Local Slope and Circulation 

From reference 4, the distributed circulation over a lifting system is 
related to the local slope by 
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where the matrix [A] Is the aerodynamic Influence coefficient matrix based on 
the paneling technique described in reference 4 . 

Circulation Specification 

Once the surface slope matrix {9z/3x} is known, chordwise integration 
can be performed to determine the local elevation surface z/c, which contains 
the effects of camber, twist, and angle of attack. The major problem to be 
solved is determining the necessary circulation matrix { T/U} to employ in 
equation (1). The problem is simplified somewhat by having the chordwise shape 
of the bound circulation remain unchanged across each span, although the chord- 
wise shape may vary from one plarform to another. The chordwise loadings allow- 
able in the program range from rectangular to right triangular toward the lead- 
ing edge and were selected because they are of known utility. An example is 
given in figure 2. Two different techniques are utilized to arrive at the span- 
wise scaling of the chordwise shapes. The particular technique to be employed 
depends on whether the configuration has dihedral. 

For a configuration having dihedral, the spanwise scaling must be deter- 
mined discretely because no finite polynomial representation of the scaling is 
known with certainty, even for ar isolated wing. However, for configurations 
with no dihedral, the spanwise scaling can be written as a polynomial for each 
planform. 


jl 1 - n , 2 ( a i + b i ", 2 + c i", 4 ) 

(see fig. 2) with a maximum of three coefficients per planform being determined 
as part of the solution. It is possible to write this polynomial as a solution 
because the isolated wing solution is known to be of the elliptical form 

'^jl - n 2, and the presence of the other planform is assumed to generate a load- 
ing disturbance which can be represented by the other two terms in addition to 
adjusting a^. Once the scaling is known from either technique, then { T/U } is 

readily obtained by multiplication. 

Lift, Pitching-Moment, and Drag Contributions 

The contributions to C T and to C , respectively, from the jth chordwise 

Li m 

row of horseshoe vortices are 



( 2 ) 
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where 


M= 4q 8 cos 6 N c 

c I*J „ I y 

m,j q S ,c , q S r c r~. 
® ref ref ref ref i«l 
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Even though 


actually occur on the wing at the jth span- 
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wise location, they can be utilized In a Trefftz plane solution If the chordwlse 
summations are performed. This utilization is possible herein because the trail- 
ing wake is assumed not to roll up, and the general configuration has specifiable 
chord loading shapes. Summing the chordwlse loadings at this point allows the 
solution of the spanwise scaling to be performed on a bent lifting line located 
in the Trefftz plane, which is, of course, ideally suited for the vortex drag 
computation. In addition, the summation reduces the number of unknowns from the 

product of N £ and to only N^. Hence, a larger value of N g can be used 

in the Trefftz plane, which should yield improved accuracy in the spanwise scal- 
ing factors without affecting the number of horseshoe vortices on the wing. 

Then, when the circulations are needed on the wing for use in equation (1), the 
well-defined variations of the spanwise scaling factors are interpolated to the 
original spanwise positions of the wing vortex lattice which is used to generate 
[A]. The procedure is implemented as follows: 

The summation in the lift expression (eq. (2)) can be written as 


N . . I . , N . . 

£S),-I «,•£ t, 

i-1 1 i-1 1 i«I+] 
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where I is the last 1 value which satisfies ■ a; that is. 



where the brackets indicate "take the greatest integer." Hence, 
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Similarly, the summation in the pitch expression (eq. (3)) can be written as 
N 
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The contribution to the vortex drag coefficient at the ith chordwise row due 
to the jth chordwise row is obtained by using only half the trailing vortex 
induced normal wash from the Trefftz plane. The result is 
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In Che ± sign, plus indicates Chat the trailing vortex filament is to the left 
of the influenced point; minus, to the right. 

In using equations (2), (3), and (10), a new vortex system is set up in the 
Trefftz plane in which the bent chord plane is represented by a system of uni- 
formly spaced trailing vortices (the quantity 2s in fig. 1). This uniformity 
of vortex spacing leads to a simplification in the equations and can be thought 
of as a discretization of the ideas of Munk (ref. 8) and Milne-Thomsen (ref. 9) 
for a bound vortex of constant strength. 

Spanwise Scaling Determination 

To determine the spanwise scaling with either technique requires the com- 
bination of the contributions from each spanwise position for configuration? 
with dihedral or the mode shape contributions for configurations without dihe- 
dral. These contributions must be employed in the appropriate total and 

C constraint equations as well as in the C n extremization operation. Due 

m UjV 


to limited space only the solution for wings without dihedral will be discussed. 
The equations to be employed in the Lagrange extremization method are 


\ C L,k 
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( 12 ) 


and 


C m " 2 ^ 6 k C m,k 
k*l 




i-1 k-1 


(13) 


where K « 6 and C. , and C , are the 

L,k m,k 

ated with the kth term in the polynomials 


°L 


and 


C contribution! associ- 
m 



(Note that k * 1, 2, and 3 are assigned to the first planform and 4, 5, and 6 

to the second.) These contributions are computed by first assuming a unit value 

of scaling with each term in the polynomial, then multiplying each resulting 

spanwise scaling distribution by the C_ . and C . terms of equations (2) 

L,j m,j 

and (3) , and finally summing spanwise over all the chordwise rows associated 
with aach set of k values (or planform). The vortex drag coefficient associ- 
ated with the ith and kth combination of spanwise scaling distributions ^ ^ 

is compared similarly. The 6^ terms are equivalent to the unknown coefficients 

in the polynomial and are the independent variables in the solution. 


An application of the preceding process to a conventional wing-tail config- 
uration is shown in figure 3. The resulting idealized loading set is of the 
type that would meet the constraints and extremlzation. 


Determination of Local Elevation Curves 


With 6, known, then {T/u}, C. , C , and C n can be determined. The 

K L 1U U|V 

results for {T/U} are interpolated to the original spanwise positions of the 
paneling which is used in equation (1) and in the following equation to find the 
local elevation curves. The equation for the local elevation above the computa- 
tional plane at a particular point (£',y) is 


f(5\?> 


f K dt 

J 1 3x 


(14) 
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RESULTS AND DISCUSSION 


General 

It is necessary to examine the sensitivity of the results of the present 
method to vortex-lattice arrangement. It is /*.lso important to compare results 
obtained with this method with those available in the literature. Unfortunately, 
the available solutions, whether exact or numerical, may not be for configura- 
tions which will exercise the constraint or extreminlzation capabilities of the 
present method. In fact, the available exact solutions are for configurations 
which are either two-dimensional sections or iscla' i three-dimensional wings 
with a nonelliptic span loading. The solutions for such configurations require 
program modifications to the span loading and involve no optimization. 


Two-Dimensional Comparison 

Various chordwise arrangements and number of vortices were investigated for 
several chordwise loading shapes, of which the a * 0.6 results are given in 
figure 4. Although difficult to see clearly from this figure, the agreement of 
the present method with analytic results (ref. 10) is good for both local slope 
and elevation. Examination of figure 4 leads to the following general conclu- 
sions_ concern '.ng the chordwise arrangement: (1) Uniform spacing is preferred; 

(2) N £ ■ 20 is a good compromise when considering both computntional require- 
ments and completely converged results. An additional conclusion is that the 
present method yields incidence angles near the leading edge which are slightly 
higher than the analytical ones. 


Number of Rows 


Along Semispan 


fs 




Various spanwlse arrangements and number of vortices were studied for one 
planform and from these studies the following conclusions were drawn: (1) Uni- 
form spacing is preferred; (2) for at least 10 spanwlse rows per semispan, the 
local slopes and elevations were not too sensitive to increasing the number. 


Three-Dimensional Comparisons 

Two comparisons with available mean-camber-surface solutions will be made. 
The comparisons are for a high-aspect-ratio sweptback and tapered wing with a 
uniform area loading at ^ ■ 1.0 and ■ 0.9V> and a lower aspect-ratio 

trapezoidal wing with a ■ 1.0, spanwlse elliptic loading at C . * 0.35, and 

M_ - 0.40. 

GO 

Figure 5 presents the predicted results from the present method for the 
sweptback wing and compares these results with those from references 1 and 11. 

A comparison of the three solutions indicates that they are all in generally 

good agreement with the exception of the results at “ 0*05. The surprising 

result is that the present method and the modified Multhopp method (ref. 1) 
agree as well as they do at this span station because of the known differences 
that exist between them near the plane of symmetry, xhe reason for the larger 


disagreement between the present method and that of reference 11 near 



0 
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is not clear, but this disagreement may be caused by the different values 

utilized by the two methods. Reference 11 effectively uses an infinite number 

since over each infinitesimal span strip across the wing the method locates a 

single quadrilateral vortex around the periphery of the enclosed area. This 

vortex extends from the leading edge to the trailing edge and includes segments 

of the edges as well. For a uniform area loading, the trailing leg parts of the 

quadrilateral vortices cancel with adjacent spanwise ones all across the wing. 

Thic leaves only the edge segments to contribute to the induced flow field. 

The present method utilizes a numerical rather than a graphical solution in 

order to provide a general capability; hence, N values ure limited as dis- 

c 

cussed previously. Also, vortices are not placed £ round the leading and trail- 
ing edges in the present method. 

A comparison of the present design method with that of reference 1 is shown 
in figure 6 for a lower aspect-ratio trapezoidal wing. The local s' opes and ele- 
vations determined by the two methods are in reasonably close agreement at the 
three spanwise locations detailed; however, an oscillatory trend is evident in 
the local slopes obtained from the method of reference 1 (fig. 6(a)). These 
oscillations apparently originate in the truncated sine series used in refer- 
ence 1 to represent a uniform chordwise distribution. Integration of the local 
slopes to obtain local elevations tends to suppress the oscillations (fig. 6(b)); 
however, the local pressures depend upon the slope rather than the elevation. 
Consequently, the measured chordwise pressure distribution will demonstrate the 
same oscillatory character. A model built according to the design of reference 1 
was tested (ref. 12), and the measured pressure distributions for a typical span- 
wise location (fig. 6(c)) indicate that indeed the oscillations are present. 
Presumably, similar measurements on a model designed by the present method would 
not behave in this manner since the input loadings are truly linear. 

Force tests (ref. 13) of an essentially identical model indicate that the 

C L 2 

measured drag polar was tangent to q that is, the vortex drag 

was indeed a minimum at the design C (or 100 percent leading-edge suction was 

L 

obtained). It is presumed from the small differences in local slope between the 
present method and the method of reference 1 that a similar result would be 
obtained for a design by the present method. 


Application to a Wing-Canard Combination 

The present method has been demonstrated by optimizing a wing-canard com- 
bination (fig. 7). To illustrate how the span load optimizing feature operates 
with the constraints, figure 8 presents individual and total span load distri- 
butions for various values of a and a., with the moment trim point at 

- c w 

Ax z c 

“0-1 a°d (This trim point is given with respect to the axis 

system shown in the sketch in figure 9.) From figure 8 there are three impor-- 
tant observations to be made: (1) The Individual span loadings change in the 

anticipated direction with the changing chord loadings in order to meet the same 
and constraints; (2) the total span loading does not change; (3) con- 

sequently, the vortex drag of the configuration is constant, as would be antic- 
ipated from Munk's stagger theorem. 
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The effects of varying the vertical separation and the moment trim point on 
the resulting drag and span loadings are also illustrated (figs. 9 and 10). All 
surfaces are designed for ^ * 0.2, a c ■ 0.6, a w * 0.8, and * 0.30 and 

have C * 0 about the moment trim point. Figure 9 shows that for all vertical 

SI 

separations, moving the moment trim point forward increases the vortex drag over 
some range, and furthermore, increasing the out-of-plane vertical separation 
reduces the vortex drag. Of course, not all moment trim points utilized will 
produce a stable configuration. These variations illustrate the importance of 
balancing the lift between the two lifting surfaces so that for some reasonable 
moment trim point and vertical separation, the vortex drag will be at a minimum. 

The minimum point on each vortex drag curve occurs with the pitching-moment con- 
straint not affecting the extremization. j 

Figure 10 presents the individual span loadings with increasing vertical ! 

separation < ® above the wing plane^ with a^ - 0.6 and a^ = 0.8. There 

are three observations which can be made from these results for increasing ver- 
tical separation: (1) The individual span loadings tend to become more ellipti- ! 

cal; (2) consequently, the vortex drag decreases; (3) the individual lift con- 
tributions show only a little sensitivity to separation distance once the canard ! 

is above the wing when compared with the coplanar results. ) 

Application to Tandem Wing Design j 

This design method has 'been employed in the determination of the local ele- J 

vation surfaces for a tandem wing. Figure 11 shows a sketch of a tandem wing 
configuration and selected results taken from the wind-tunnel tests made with a f 

model based on this design at a Mach number of 0 30 (ref. 14). At s 0.35 • 

the vortex drag increment is correctly estimated. The measured C is slightly 

® | 

positive (0.02). Reference 14 states that a part of the error ^C m should be 

zero) is a result of a difference in the fuselage length between the designed and 
constructed model. ! 


Design of a Wing-Winglet Configuration 

Figure 12 presents the wing-winglet combination of interest along with 
pertinent aerodynamic characteristics and local elevations obtained from the 
present method. For comparison these same items are calculated with a program 
modification that adds a root-bending-moment constraint to produce the same 
moment that would be obtained on the original wing extending to the plane of 
symmetry but without its basic wingtip. The assumed span loading is elliptical. 
The force and moment coefficients are based on the wing outside of a representa- 
tive fuselage and without the basic wingtip. 

The results of this comparison are as follows: (1) The root-bending-moment 

constraint increases the vortex drag slightly because of the changes in the c^c 

distribution required; (2) the differences in local elevations are confined pri- 
marily to the outer 50 percent semispan and result mainly from the differences 
in the incidence angles; (3) significant amounts of incidence are required in 
the winglet region with or without the root-bending-moment constraint. 
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CONCLUDING REMARKS 


A new subsonic method has been developed by which the mean camber (local 
elevation) surface can be determined for trimmed noncoplanar planforms with 
minimum vortex drag. This method employs a vortex lattice and overcomes pre- 
vious difficulties with chord loading specification. This method designs con- 
figurations to have their local midsurface elevations determined to yield the 
span load for minimum vortex drag while simultaneously controlling the pitching- 
moment or root-bending-moment constraint at the design lift coefficient. This 
method can be used for planforms which (1) are isolated, (2) are in pairs, 

(3) include a winglet, or (4) employ variable sweep, but only at a specified 
sweep position. 

Results obtained with this method are comparable with those from other 
methods for appropriate planforms. The versatility of the present method has 
been demonstrated by application to (1) isolated wings, (2) wing-canard config- 
urations, (3) a tandem wing, and (4) a wing-winglet configuration. 
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Figure 3.- Idealized loading set on trimmed configuration for minimum drag. 
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Figure 6.- Local slopes, elevations, and lifting pressure distribution! 

C, . - 0.35; M_ - 0.40. 



Figure 7.- Typical wing-canard combination. 
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Figure 8.- Effect of chord loading on span loadings for trimmed 
coplanar wing-canard c~mbination; M ■ 0.30. 
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Figure 9.- Vortex drag for range of center-of-gravity positions and vertical 

separations; C T , * 0.2; M_ * 0.30. 
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Figure 10.- Effect of vertical displacement of span loadings for trimmed wing- 


canard combination 
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Figure .11.- Longitudinal aerodynamic characteristics of tandem wing; 

C. . - 0.35; M - 0.30. 
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(a) Aerodynamic characteristics. 





(b) Incidence angle distribution. 
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(c) Local elevations. 


Figure 12.- Effect of root-bending-moment constraint on aerodynamic charac 
teristics, local elevations, and incidence angle distribution of wing- 
vinglet combination (AR ■ 6.67); N ■ 20; N “17; M„ * 0.80. 



